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We review various cosmological models with a local underdense region (local void) and the aver- 
aged models with the backreaction of inhomogeneities, which have been proposed to explain (with- 
out assuming a positive cosmological constant) the observed accelerating behaviors appearing in 
the magnitude-redshift relation of SNIa. To clarify their reality, we consider their consistency with 
the other observational studies such as CMB temperature anisotropy, baryon acoustic oscillation, 
kinematic Sunyaev-Zeldovich effect, and so on. ft is found as a result that many inhomogeneous 

■ models seem to be ruled out and only models with the parameters in the narrow range remain to 

be examined, and that, unless we assume very high amplitudes of perturbations or gravitational 

\ energies, the averaged models cannot have the accelerated expansion and the fitted effective A has 

not the value necessary for the observed acceleration. 

PACS numbers: 98.80.-k, 98.70.Vc, 04.25.Nx 
; I. INTRODUCTION 

' The observational magnitude (m) - redshift (z) relation for distant supernovas of type la (SNIa) have been 
studied by the High-z SN Search Team[ll, 0, 0] and the Supernova Cosmology Project [^], and among the 

[ Friedman-Lemaitre-Robertson- Walker cosmological models, the A-dominated models have been selected as the 
best models which can explain their accelerating behavior. At present, the model parameters (SImj^^a) — 
(0.25, 0.75) are regarded as the representative one. 

About 40 years ago, the cosmological constant havs been paid attention to avoid the contradiction between 
the age of oldest stars and the cosmic age[l,[l,0,[l], but its value has been uncertain, since the Hubble constant 
was inaccurate. But recently the existence of A has been realistic and reduced to the above value. 

The A-dominated models with the above parameters have recently been found to be consistent with all 
results of many observations such as CMB temperature anisotropies[§, [lol |. the baryon acoustic oscillation 

'■ (BAO)[ll[il,lll,[il[il, statistics of cluster abundance [H, [ll and the inverse Sachs-Wolfe effect [H, HO] , 
and so it is regarded as a concordant model. 

In these models, however, there are well-known unsolved problems: the cosmological-constant problem and 
the incidence problem. The former problem shows that the standard value of the cosmological constant is 
too small (~ 10~^^°), compared with a constant appearing as the zero-pont fluctuation of a field. The latter 
shows that the role of the cosmological constant is important only at the later stage of z < 1, when structure 
formations are remarkable. Moreover we know the low-Z anomaly in the CMB temperature anisotropics which 
remain to be unsolved [2ll. [2^ . 

In view of these problems, on the other hand, astrophysical explanations for the accelerating property of SNIa 
were tried by considering the optical role of inhomogeneous matter distribution which may cause the apparent 
cosmological acceleration without A or dark energy. The explanations are divided into the following two types 
by what inhomogeneity is assumed: 

' (1) Non-Copernican inhomogeneity. We assume to live in a special point of the zero A universe with a spherically 
symmetric underdense local inhomogeneity, which is called a local void. The apparent accelerating behavior of 
SNIa is obtained by deriving the m-z relation in this inhomogeneous model. 

(2) Copernican inhomogeneity. The universe is assumed to have an uniform distribution of density perturbations 
with zero A. The accelerating behavior is obtained by deriving the effective A from them in the two procedures: 

a) Averaging and backreaction of inhomogeneous models with A = 0. First the average model with the effective 
A is derived and the m-z relation is obtained to to examine the accelerating behavior. 

b) Fitting. By comparing the m-z relation in a perturbed model with zero A and that in an unperturbed model 
with nonzero A, the best fitted value of effective A is obtained. 

The excelent review about the works on (1) and (2)a), b) was written by Celerier[23j in 2007. In this note 
we have a review of the works which include more recent theoretical and observational contributions. 
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II. A LOCAL INHOMOGENEITY AND SNIA 

Soon after the contributions by the High-z SN Search Team[li[li3 and the Supernova Cosm olog y Project Hi , 
another possibihty was suggested independently by Celerier[2J|, Goodwin et al. (unpublished) [25| . and nie|26l. 
[27} . It was that the m-z relation of SNIa may be reproduced also in non-Copernican inhomogeneous cosmological 
models with a central underdense region. 

Celerier[2^ discussed the accelerated expansion from the qualitative viewpoint with general inhomogeneous 
models, and Goodwin et al.|25j] discussed at the similar epoch with a physical analysis of SNIa result due to the 
local to global Hubble-constant ratio. 

In my first model(26l. |27|. the inhomogcncity is described by an open FRW solution (in the inner underdense 
region) and a self-similar solution (in the outer overdense region). The latter is a special case of Lemaitre- 
Tolman-Bondi (LTB) solutions which is connected with the former open FRW solution at a boundary, in order 
to avoid the central singularity, and tends to the Einstein-de Sitter (EdS) solution in the limit of radial infinity. 
It was found in this model that the theoretical m-z relation can be consistent with the observed relation, by 
adjusting the Hubble constants (Hq) and the total density parameters {i^Mo) in the two regions, in the case 
when the boundary is at the distance of z ~ 0.07 from the center. The super-horizon version of this model, in 
which the distance of the boundary is larger than the horizon size, was once proposed in 1995[28j. 

In my second model[2^, two FRW models with A = are used to build an inhomogeneous model, that is, 
an open FRW model (in the inner underdense region) and a flat FRW model (in the outer overdense region) 
connected with a discontinuous, sharp boundary with z « 0.07. The angular diameter distance dA is derived, 
solving the Dyer-Roeder equation, the theoretical m-z relation was derived, and it was shown that the observed 
m-z relation can be reproduced similarly by adjusting the values of Hq and ^mo in the two regions. 

Moreover, in my next paper [30|. the quantitative comparison between the observed and theoretical m-z 
relations was done using the second inhomogeneous model and the best model parameters were estimated. 

Subsequently, Iguchi, Nakamura and Nakao 331] analyzed general inhomogeneous models using LTB solutions. 
They reproduced the m-z relation in the concordant model, and it was shown that, at an intermediate radius 
(with z < 1.7), a critical point appears in the models when we assume the uniform big-bang time and asymptotic 
vanishing spatial curvature. 

The geometrical structure of LTB models mimicking A or dark energy was analyzed by Vanderveld et al. [33 | 
in detail and it was shown that a weak central singularity and a critical point appear generally. Recently Yoo, 
Kai and Nakao[33| have showed that the LTB models, which reproduce the m-z relation in the concordant model 
and have no critical point, can be obtained only under the condition of uniform big-bang time, by improving 
Iguchi et al.'s model, and Clifton et al. 34j| derived LTB models without central weak singularity and critical 
point and discussed the reproduction of the observed m-z relation in their model. 

The works on reproducing the observed m-z relation in LTB models were done successively by Alnes, 
Amarzguioui and Gr0n 35, 36|, Alnes and Amarzguioui.33, [1^, MansourifH, H^], Moffat [ill, li^ . Biswas et 
al.jill, and Alexander et al.|44]. 

On the other hand, Kasai|i45j| found from SN data themselves @that they can be divided into the low z group 
(with z < 0.2) and the high z group (with z > 0.3), which correspond to higher and lower Hubble constants, and 
that the different trend of the data with respect to redshifts may represent the inhomogeneity of cosmological 
models. 



III. CONSISTENCY OF INHOMOGENEOUS LOCAL- VOID MODELS WITH THE OTHER 

OBSERVATIONS 

In order that the inhomogeneous local- void cosmological models may be realistic, they must be consistent with 
not only SNIa data, but also the other observations such as CMB temperature anisotropics, BAO, the kinematic 
Sunyaev-Zeldovich effect, and so on. In the following, let us review the works investigating the consistency with 
these observations. 



A. CMB temperature anisotropics 

We can see several acoustic peaks in the correlation (Cj) - multipole (l) diagram for I > 200, given with the 
WMAP data. Since the inhomogeneous models are equal or nearly equal to the Einstein-de Sitter (EdS) model 
at the recombination epoch, the EdS model must be consistent with the observed property of CMB anisotropics. 
First, Alnes et al.[35|] discussed the consistency for the first peak and Alexander et al.Q showed that their 
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model (the Minimum model) could give the consistent first and second peaks by assuming a value of the Hubble 
constant (in the outer region) significantly lower than the conventionally accented value (70). The consistency 
with CMB was discussed also by Blanchard et al.46] and Hund and Sarkar|43 in small values of the Hubble 
constant (~ 47) in the outer region. 

If an observer is off-center, he observes the dipole component of CMB temperature anisotropics, proportional 
to the distance r from the center of the models. So this value r is constrained by the observed upper limit of 
the dipole anisotropy and the upper limit of r is about 15 Mpc [ssl. l4ll. Ilsj . 

B. BAO 

The above acoustic peaks (representing an oscillation of photon-baryon fluid around and before the recombi- 
nation epoch) appear also as the baryon acoustic oscillation (BAO), which is an oscillation of baryons (at the 
later stage) imprinted in the matter spectrum. 

Recently it was detected in the SDSS and 2dFGRS surveys, and the BAO scale has been regarded as a 
standard ruler pH [T^ . [Tsl . [Tsj . The characteristic (BAO) scale is the sound horizon at the recombination 
epoch : 

/"OO 

TsiZrec) = dz Cs{z)/H(z), (3.1) 

J Zrec 

where Cs{z) is the sound speed at redshift z and z^ec is the redshift at the recombination epoch. This scale is 
approximately expressed as 

r, w 147(f7o/iV0-13)~°-^(r2b/iV0-024)-" "* Mpc, (3.2) 

where the density and Hubble parameters have the values at the place we notice. 

The observed scales of BAO brought from the galaxy samples in the above surveys are used to constrain the 
values of the distance measure dy(z) = [(1 -I- z)^(d^)^cz/if(z)]-^/'^. Here dA is the angular diameter distance 
and H{z) is the Hubble parameter at redshift z, and dv is the average of the scales perpendicular and parallel 
to the line-of-sight directions. Percival et al.[l^ have recently obtained the following relations at z = 0.2 and 
0.35: 

rs/dv{0.2) = 0.1980 ±0.0058 
rs/dv{0.35) = 0.1094 ±0.0033. (3.3) 

In the A dominated concordant models, rs/dy(0.2) and rs/dv{0-35) can be reproduced approximately (in about 
95%). In my second models with a local void, on the other hand, these two ratios are about 80% of the observed 
values. This is because the epochs z = 0.2 and 0.35 belong to the outer region with the EdS model. So my 
second models cannot reproduce them, and they are therefore ruled out. In my first model and the Minimum 
model given by Alexander et al.[4^ also, the situation is similar and they are ruled out, as long as the boundary 
is of the order of 300 Mpc. 

In order that inhomogeneous models may not contradict with the BAO observation, the epoch z = 0.35 must 
belong to the inner underdense region, so that the scale of the inhomogeneous models must be Gpc size. From 
this viewpoint, the Gpc-size inhomogeneous (LTB) models have recently been studied by Clifton et al.fU and 
Garci'a-Bellido and Haugb0llef49|. 

It has recently been found that the scale of BAO in the parallel to line-of-sight or in the radial direction 
(the radial BAO scale) impose on the models a more stringent condition than the BAO in the perpendicular 
direction. RBAO has been studied by Gaztahaga et al.flO] and it was found that the concordant models are 
consistent with the observational result of RBAO. On the other hand, Zibin et aH51j] showed that RBAO 
imposes stringent conditions to Gpc-size inhomogeneous models, using Gaztanaga etal.'s data at z = 0.24 and 
0.43. They considered the two types of models: a constrained model and an unconstrained model, which were 
made so as to reproduce the recent SN data and WMAP data with the boundary of z « 1, and investigated the 
consistency with RBAO for the above two redshifts. The constraint condition is / Sp{ti)r'^dr < 0. Their result 
is 

(1) the unconstrained model is consistent with the RBAO data, but Hq (the present Hubble constant) must 
be ~ 44, so that this model is ruled out, and 

(2) the constrained model is approximately consistent with the RBAO data and Hq is ^ 60, so that this model 
is not ruled out at present. However the consistency is not so good as the concordant models. 
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C. Kinematic Sunyaev-Zeldovich effect 



The Sunyaev-Zeldovich (SZ) effect is a small spectral distortion of the CMB radiation caused by the collisions 
of CMB photons with hot thermal electrons. CMB photons passing through the hot center of massive clusters 
interact with their electrons and take a small distortion in the CMB spectrum due to the inverse- Compton 
scattering. 

If the clusters are moving with respect to the CMB rest frame, there is an additional spectral distortion (the 
kinematic Sunyaev-Zeldovich effect) due to the Doppler effect of the cluster velocities on the scattered CMB 
photons. If Vpec is the component of the cluster velocity along the line of sight, then the Doppler effect leads to 
the following distortion of CMB spectrum: 

^TeiVpec/c), (3.4) 

where Te = neCrR {ne, o't, R are the electron density, the Thompson scattering cross-section and the effective 
radius of a cluster) . 

In inhomogeneous models with a local void, the CMB photon received by a central observer is emitted at 
the recombination epoch in the outer EdS region, in which the Hubble constant (iJo)eds is ~ 47 km/s/Mpc. In 
the regions inside and in the neighbourhood of the local void, the Hubble constant {Hq)ioc is assumed to be 
larger than {Ho)eds- So a cluster in the distance r from the center of the inhomogeneous models has a velocity 
[{Ho)ioc — {Ho)eds]r, relative to the CMB rest frame. This velocity can be observed as a peculiar velocity (wpec) 
of clusters in the kinematic Sunyaev-Zeldovich effect. 

The possibility to observe the cluster velocities systematically and put a constraint on cosmological models 
was studied by Benson et al.[5l|- Recently, Garcfa-Bellido and Haugb0llej53] have shown the constraint due to 
the kinematic Sunyaev-Zeldovich effect using observed 9 clusters, and shown that a strong constraint is given to 
the Gpc LTB models. It is found that the models with the local void region larger than ~ 1.5 Gpc are ruled out 
practically, and only special Gpc LTB models with a limited range of cosmological parameters may be allowed. 
But, if a systematic peculiar motion of clusters is discovered, this effect may support the local-void model. 



D. Spectral distortion of CMB radiation in the reionized region 



In inhomogeneous models with a local void, the inner void region is considered to be at the reionized stage. 
The ionized gas there is moving outward, relative to the CMB frame, and leads to the Doppler effect. When 
CMB photons are reflected by this ionized gas and reach an observer at the center of the void region, spectral 
distortion appears in the accepted CMB radiation owing to the Doppler anisotropy. Here the reionized matter 
plays a role of moving mirror. Caldwell and Stebbins 54'| derived this distortion in the inhomogeneous models 
with a local void of various sizes and showed that the constraints for the models can be obtained from the 
observed upper limit to the distortion of CMB spectrum. As a result, they found that severe constraints can 
be obtained and the models with the local void of the largest size (~ 2.5 Gpc) are ruled out. 



IV. UNIFORM DISTRIBUTIONS OF DENSITY PERTURBATIONS 



The appearance of observational acceleration due to the Copernican uniform distribution of density pertur- 
bations has been studied by the following two different forms. 



A. Averaging and backreaction of inhomogeneous models 



include Buchert[55j, Buchert 
M [el and Nambulei 



and 
M 



Averaging and backreaction have been studied by many workers, who 
CorforajHl, Ellis and Buchert^], Kolb et alpl, [H, [13] , Kolb et al.|6l|, Kasai[6 
[13, [fill for various gauges, and Zalaletdinov[69j| and ParanjapejzOl for a covariant form. 

In the Buchert formalism, we consider an inhomogeneous model with dust, and its metric is expressed in the 
comoving and synchronous gauge as 



-dt^ + q,j{t,x"^)dx'dx^. 



(4.1) 
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For averaging, we specify a hypersurface S of constant t, and take a compact region D of E. If the volume of 
D is Vd and -0 is a scalar variable, the average < tp >]:, of ip over D is 

< ^ >D= / ^dS. (4.2) 

yo Jd 

Especially, for the matter density we have 



<P>D=^ I pdY.. (4.3) 



The averaged scale factor is defined by 



ao = {VdY'\ (4.4) 



Then we obtain from the Einstein equation 



— ^ < P>D +Qd, 

< p>D <n>D -^Qd (4.5) 
and 

ia%QDy +a%ial<n>Dy =0, (4.6) 

where TZ is the scalar curvature of E, 

Qd = ^(< e^>D-<0 >l) - (aya'J)i3, (4.7) 

and 9 is the expansion of the world line of the dust fluid. In the derivation of these equations, the incommutability 
between the time derivative and the averaging was used: 

< V >D=< tp >D + < >D - < >D< i/j >D, (4.8) 

where < ip >d=< dip/dt >o and < ip >jj= d < ij] > d /dt. The condition of averaged acceleration du > is 
given by 

Qd >^k^ <P>d ■ (4.9) 

If this condition is satisfied, it seems that the averaged acceleration may be realized, and so it has been studied 
by many workers under what situation it is satisfied. 

Nambu and TanimotofTll considered the region D which consists of many homogeneous and isotropic small 
regions with different scale factors. Then it was shown that, even if we have deceleration in each region, the 
average scale factor ajj can have acceleration such as du > 0. That is, the deceleration observed in each region 
can be compatible with the averaged acceleration. Ishibashi and Wald[7^ discussed this situation and found 
that the Buchert averaging procedure has ambiguity both with regard to the choice of time slicing and the 
choice of domain Z?, and that their special choices can artificially derive the averaged acceleration. 

Independently of Nambu and Tanimoto, Rasanenj?^ considered a model with two disjoint regions consisting 
of an overdense region and a completely empty region. In both regions we have separately FRW spacetimes 
with scale factors ai and 02, where ai oc (1 — cosm) and t cc (u — sinw) and 02 oc t. If we define the averaged 
scale factor a by a'^ = ai^ + 02^, then the averaged deceleration parameter q = (da)/d^ becomes negative 
and so acclerating, when the overdense region turns around and starts collapse. However, Rasanen's model is 
physically incomplete, because the junction condition is neither used nor satisfied. Paranjape and Singh[8l|. 
on the other hand, derived a LTB solution corresponding to Rasanen's model, in which there are three regions 
corresponding to Rasanen's two regions and a medium region. As a result, they found that q < cannot be 
realized and so we have not obtained the averaged acceleration in the physically adjusted solution. 

Next, Ishibashi and Wald[73| insisted that in the Newtonian gauge the perturbed universes can be expressed 
by the metric 



Van- 



ds^ = -(1 + 2'i>)dt^ + a^{t){l - 2^)j,jdx'dx^ , 



(4.10) 
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where 7^ denotes the metric of constant-curvature space. The potential 4" is related to the matter density by 
the cosmological Poisson equation 

a-2A(3)Vl/=i«;2^p=i«2(^_^)^ (4.11) 

where A(3) = and Di denotes the derivative corresponding to 'jij. Here '3/ satifies \'^\ << 

1, << a-^iP^^'A^', and << {D' 0^)0,0^"^ . In this expression of spacetimes, the devia- 

tion of perturbed universes from the FLRW spacetime is very small, even if there are nonlinear perturbations 
like galaxies, clusters, voids and superclusters. This is because ^ ~ G6M/{c^R) << 1 for such perturbations, 
where 6M « 6pR^ and R is the radii of these structures. 

The above analyses in the Newtonian gauge have a long history. Nariai and Ueno[74] and Irvine 75] derived 
cosmological equations in the cosmological Newtonian approximation, and they have been applied to the nonlin- 
ear treatment of matter evolution and the N-body simulation. Subsequently the formulation in the cosmological 
post-Newtonian approximation also was treated|76l FttI pTOj . 

From the viewpoint of Newtonian approximation, it is expected that the averaging and backreaction of sub- 
horizon perturbations have only very small contributions to the background. Kasai et aH80| derived the no-go 
theorem in the Newtonian gauge that the nonlinear backreaction neither accelerates nor decelerates the cosmic 
expansion, because the cosmic averaged acceleration d/a is determined merely by the mean density as 

a 47rG' 

- = ^<P> (4.12) 

and the nonlinear backreaction reduces the expansion rate d/a as 

(^)' = ^<P>-^<i^W*> (4.13) 

Paranjape and Singh[8l| and Siegel and Fry|82,] also analyzed the possibility of averaged expansion in the 
Newtonian gauge and obtained the negative results. 

In the super-horizon case, Kolb et al.fH, [gOl and Barausse et al.fsSl studied the possibility of averaged 
acceleration which is caused by second-order super-horizon perturbations associated with primordial inflation, 
but Flanaganjs^l and Geshnizjani et al.jsl] obtained the negative results because of the incompleteness in their 
second-order analysis and the constraints from the CMB anisotropics. 

On the other hand, Kai et al.[86] studied the possibility of accelerated expansion by use of the LTB solution 
and found that, in a strong inhomogeneous condition such as in Rasanen's one, the acceleration of cosmic volume 
expansion is realized only in some cases when the size of the perturbed region is comparable to the horizon 
radius of the EdS universe. 

From the above studies, it seems difficult, therefore, that the accelerated expansion is caused not only by the 
averaging and backreaction of sub-horizon perturbations, but also by those of super-horizon perturbations. 

Wiltshire [13, [H, [s^ studied the difference of gravitational energies and their influence to the clock rates in 
the different regions such as the averagedly expanding region, voids and walls (finite infinity), and considers 
that this difference brings the observed accelerated expansion. Here we take notice of the definition of the 
gravitational energy and why it can be so large as to change remarkably the clock rates in different regions. 



B. Fitting 

Without averaging and backreaction of inhomogeneous models, we can consider the accelerating effect of 
inhomogcncities by calculating directly the observational quantities such as the angular diameter distance and 
the redshift in inhomogeneous models and comparing them with the counterparts in FLRW models with nonzero 
A. That is, by fitting these two models, we can obtain an effective A. 

Vanderveld et al. 9(i| studied this fitting by using the post-Newtonian approximation and assuming the linear 
perturbations which are normalized with respect to CMB temperature anisotropics and cluster statistics. Then 
they obtained an important result that the effective A is about 0.004, which is very small, compared with the 
value required by the realization of the observed acceleration. 

Recently, Marra et al.[91, 92] studied the fitting by assuming nonlinear perturbations whose amplitudes are 
much larger than those expected from the CMB normalization and which are given using many arranged swiss- 
cheese models on scales of several hundred Mpc. From these models they could obtain the effective A which 
is comparable with the value necessary for the observed acceleration. At present, however, we do not know 
how these perturbations with such high amplitudes can exist, without contradicting with the observed CMB 
temperature anisopropies and the other obsevations. 
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V. CONCLUDING REMARKS 

Various inhomogeneous cosmological models with a local void have so far been proposed to explain the ob- 
served accelerating property in the m-z relation of SNIa, without a cosmological constant, and their geometrical 
structure of the models obtained using LTB solutions also has been studied in details. At present, on the other 
hand, there is no data for SNIa with z > 1.7. We should notice that this is very important for the model 
selection, especially to clarify which of the concordant model or local-void models are better with respect to 
this observation. 

Recently the studies of their consistency with observations such as BAO (especially RE AS), kinematic 
Sunyaev-Zeldovich effect, and spectral distortion at the reionization stage have shown that many models with 
a local void are ruled out, and we find that only Gps models with the narrow range of parameters remain to be 
examined. So at present the possibility for these models to survive is small, but may not be zero. 

Moreover, the averaging and backreaction of inhomogeneous models and the fitting with nonzero A models 
have been studied. At present, however, it seems difficult to obtain accelerated expansion and the expected 
effective A, unless we assume perturbations with amplitudes much larger than those corresponding to the CMB 
normalization or high amplitudes of gravitational energies included within structures like clusters and voids. 

In these analyses, we have assumed the models based on the Einstein gravitational theory and the existence 
of an inflationary early stage. If the models should be derived, however, from the other gravitational theories, 
such as the superstring theory, the cosmological situation may be quite different, because we do not know how 
the inflation arises, what inflation we have, and whether the expected cosmological constant can exist or not. 
Then the inhomogeneous models with a local void may play some more role to explain the observed accelerating 
behavior. 
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